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We demonstrate that the addition of a single carbon impurity leads to significant 
changes in the thermodynamic properties of Ni clusters consisting of more than 
a hundred atoms. The magnitude of the change induced is dependent upon the 
parameters of the Ni-C interaction. Hence, thermodynamic properties of Ni clusters 
can be effectively tuned by the addition of an impurity of a particular type. We also 
show that the presence of a carbon impurity considerably changes the mobility and 
diffusion of atoms in the Ni cluster at temperatures close to its melting point. The 
calculated diffusion coefficients of the carbon impurity in the Ni cluster can be used 
for a reliable estimate of the growth rate of carbon nanotubes. 



I. INTRODUCTION 

Thermodynamic properties of atomic clusters depend considerably on the cluster size. It 
has been shown that the melting temperature of a small spherical particle decreases with the 
reduction of its radius .-'^'^'^'^ This is due to the substantial increase in the relative number 
of weakly bounded atoms on the surface in comparison with those in the bulk. Such a size 
effect on the melting temperature of small particles of different metals having diameters 
down to 2-3 nm has been confirmed in series of experiments.-'^*^'^ 

However for clusters having sizes smaller than 1-2 nm, the melting temperature is no 
longer a monotonic function of the cluster size. Experiments on sodium clusters NaTv, with 
number of atoms = 50 — 360, have demonstrated that melting temperature as a function of 
size shows a prominent irregular structure with local maxima.— >iiii^ii^ The origin of the non- 
monotonic variation in the melting temperature with respect to cluster size lies in interplay 
between electronic and geometric shell effects in the sodium clusters.— Intensive theoretical 
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efforts have been undertaken to identify the details of the geometric and electronic structures 
underlying the variations in the melting temperature.— 

Experiments on small cluster ions of tin^^ and galliunt^^ have confirmed the violation of 
the linear relationship between the reduction in the melting temperature and the inverse 
radius of the cluster. It was discovered that the melting temperature of selected Sn^v and 
Ga^v clusters, of sizes < 40, can considerably exceed the melting temperature of the 
corresponding bulk material.— This behavior was explained by the structural differences 
between the small clusters and the bulk.— 

Modification of a cluster structure can also be induced by the doping of an impurity or 
substitution of one or several atoms in a pure cluster with an atom of a different type.— 
It has been shown that doping of sodium clusters with Li, Cs and O2 impurities (as well 
as AI13 and Gai3 clusters with a carbon impurity) results in a decrease of the melting 
temperature of the cluster.— Whereas doping of the icosahedral silver clusters with Ni 
and Cu atoms considerably increases the melting temperature.— Therefore one can suppose 
that the thermodynamic properties of a selected cluster can be tuned by doping with an 
impurity of a particular type. 

In this paper we report the results of a systematic theoretical study regarding the effect 
of impurity on the thermodynamic properties of Ni clusters. We demonstrate that adding a 
single carbon impurity can result in changes in the melting temperature of an Nii47 cluster. 
The magnitude of the change induced is dependent upon the parameters of the interaction 
between the Ni atoms and the C impurity. Hence, thermodynamic properties of Ni clusters 
can be effectively tuned by the addition of an impurity. We also show that the presence of a 
carbon impurity considerably changes the mobility of atoms in the Ni cluster at temperatures 
close to its melting point. 

The choice of Ni clusters is stipulated by their high chemical and catalytic reactivity, 
unique properties and multiple applications in nanostructured materials.— An important 
example of such an apphcation is the process of the catalytically activated growth of carbon 



nanotubes. The mechanism of this process is not yet well understood (see, e.g., Refs. 



31 



30, 



32I and references therein) and knowledge of the specific role of the impurity in the Ni 



catalytic nanoparticle may ascertain whether the carbon nanotube structure and its growth 
kinetics can be controlled. The thermodynamic state of the catalytic nanoparticle plays a 
crucial role on the carbon nanotube growth. '^^ The nanotube growth rate can be obtained 



3 



by a solution of a set of kinetic equations which include, in particular, the diffusion flux of 
carbon through the metal particle.— On the other hand the diffusion coefficients depend 
on the thermodynamic state of the catalytic Ni cluster which can, in turn, be tuned by doping 
with an impurity. Therefore direct molecular dynamics study of carbon diffusion in the Ni 
cluster is important for a reliable estimation of the growth rate of carbon nanotubes. 



II. THEORETICAL METHODS 



The study of structural and dynamical properties of the clusters of transition metals is a 
challenging task due to the presence of unfilled valence (i-orbitals. The high density of the 
c?-states and their localized character make the direct ab initio methods computationally 
very demanding for clusters larger than several dozens of atoms.— In order to describe 
the structure of clusters of larger sizes, one needs to use approximate methods and model 
interatomic potentials. 

An effective approach for study of transition-metal clusters is the embedded-atom 
method^Ii^i^iiSdid^ which takes into account many-body effects. The latter appears through 
the inhomogeneous electron density of the system. In this paper, the molecular dynamics 
study of the Nii47 cluster has been performed using the Sutton-Chen^ many-body potential 
which belongs to the family of the embedded-atom types of potentials. The Sutton-Chen 
potential has been shown to reproduce bulk and surfac e prop erties of transition metals and 
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46l and references therein). 



their alloys with sufficient accuracy (see, e.g., Refs. 
The potential energy of the finite system within the Sutton-Chen model has the following 
form: 
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Here rij is the distance between atoms i and j, e is a parameter with dimension of energy, 
a is the lattice constant, c is a dimensionless parameter, and n and m are positive integers 
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with n > m. The parameters provided by Sutton and Chen for nickel have the following 
values:.^ e = 1.5707 ■ lO'^ eV, a = 3.52 A c = 39.432, n = 9, and m = 6. 

The determination of a reliable model potential for the Ni - C interaction is a difficult 



task. In Refs. 
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48 



49 



the binding energy and different charge state of various forms of NiCn 
clusters (n= 1-3) were determined using ab initio density functional calculations based on 
the three-parameter Becke-type gradient-corrected exchange functional with the gradient- 
corrected correlation functional of Lee, Yang, and Parr (B3LYP)^i^i^ and LANL2DZ basis 
set (see, e.g., Ref. 5^ and references therein). This was done in order to develop a many- 
body potential for the transition metal - carbon interaction. The many-body potential was 
then successfully used for a molecular dynamics study of the formation of metallofuUerenes^ 
and single-walled carbon nanotubes.— 

In the present work, however, the interaction between Ni atoms and C impurity is modeled 
by the Morse pair potential: 

V^'-^{r) = Em ((1 - e''(i-'-/'-°))' - l) . (3) 

We have determined the parameters of the Morse potential [em = 2.431 e V, p = 3.295, 
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within the 



tq = 1.763 A ) by fitting the Ni - C interaction obtained in Refs. 
B3LYP/LANL2DZ method. The choice of using the pairwise Morse potential is due to its 
simplicity which allows us to study the infiuence of the parameters in ([3]) on the thermo- 
dynamic properties of the C-doped Nii47 clusters, keeping a clear physical picture of the 
process occuring in the system. 

The optimized initial geometries of the cluster have been determined by finding local 
minima on the multidimensional potential energy surface. We have applied an efficient 
scheme of global optimization, called the Cluster Fusion Algorithm (CFA).— i^"^ The CPA 
belongs to the class of genetic global optimization methods,—"^ which are very promising for 
structural optimization of nanoalloys, see, e.g., Ref. |59| and references therein. The scheme 



has been designed within the context of determination of the most stable cluster geometries 
and it is applicable for various types of clusters.— i^*^*^ 

Molecular dynamics simulations have been performed for the canonical (NVT) ensemble 
of particles: the number of particles N, the volume V and the temperature T of the system 
are kept constant. Integration of Newton's equations of motion have been done using the 
Verlet leapfrog algorithm, with a time step At = 1 fs and a total simulation run of 10 ns 
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(excluding an initial equilibration time of 50 ps). For temperature control, we have used the 
Nose-Hoover thermostat^i^ because it correctly generates a canonical ensemble. 

III. NUMERICAL RESULTS AND DISCUSSION 



During the last, decade numerous theoretical and experimental efforts have been devoted 
to studying structural, electronic, and magnetic properties of Ni clusters (recent compre- 



hensive reviews can be found in Refs. 



65 



66l ). It has been demonstrated that for small Ni 



clusters, a motif based on the icosahedral structure dominates the cluster growth, at least 
around sizes of complete Mackay icosahedra.— '^'^'^ The results of the geometry optimiza- 
tion for the pure and the C-doped Nii47 clusters are shown in Fig. [1] The ground state 
of the Nii47 cluster is a perfect icosahedron. Fig. [It^a), in agreement with the results of the 
previous works (see Refs. I65lj69l and references therein). 




a) 



b) 



c) 



FIG. 1: (Color online) — (a) Optimized structure of a pure Nii47 cluster; (b) the isomer state 
structure of the C-doped Nii47 cluster; (c) the ground state structure of the C-doped Nii47 cluster. 



For the C-doped Nii47 cluster we found that the lowest energy state possesses the icosa- 
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hedral type of structure with the C atom in the central position and one Ni atom outside the 
filled icosahedral shell, as shown in Fig.[T](c). The energetically closest isomer of the C-doped 
Nii47 cluster is the deformed icosahedron with the C impurity located in the vicinity of the 
cluster center, as depicted in Fig. [TJ^b). The impurity creates a local distortion between the 
first and the second closed icosahedral shells. The difference in the binding energy per atom 
between the ground. Fig. [U^c), and the isomer. Fig. [T](b), states of the C-doped Nii47 cluster 
is 0.002 eV. The isomer state. Fig. [TJ^b), is metastable because in order to transform the 
cluster structure from the isomer. Fig. [D^b), to its ground state. Fig. [U^c), it is necessary to 
overcome the energy barrier. The isomer. Fig. [I](b), can be naturally formed by successive 
migration of the outer C atom from the cluster surface towards its center. Namely, this type 
of isomer are presumably formed in carbon nanotube growth experiments where a molecule 
of a feedstock gas or a free C atom collides with the surface of the Ni particle. Therefore, 
we have performed molecular dynamics simulations for the initial geometry corresponding 
to the isomer state. Fig. [U^b) of the C-doped Nii47 cluster. 



A. Pure Nii47 clusters 

Figures [2] and [3] demonstrate the temperature dependence of the time-averaged total 
energy ( Etot ) (i-e. the caloric curve) and the heat capacity Cy calculated for the Ni^y 
cluster, respectively. The heat capacity at constant volume is defined as a derivative of the 
internal energy over the temperature: Cy = {dE/dT)y. 

A thermal phase transition is indicated in the caloric curve by a change in the gradient 
of the temperature-dependent total energy (Etot)- The height of the jump near the phase 
transition point gives an estimate of the latent heat, i.e. the energy which is associated 
with the destruction of the ordered lattice. Figure [2] demonstrates that the change in the 
slope of the caloric curve takes place in the wide interval of temperatures between T freezing 
~ 600 K and Tmeiung ~ 800 K. For the temperatures T < T freezing, the cluster is completely 
frozen; while for temperatures T > Tmeiung, the cluster has melted. The intermediate interval 
T freezing < T < Tmeiting correspouds to the mixed state where the solid and the liquid phases 
coexist. This is a typical behavior of the caloric curve for finite systems, in contrast to the 
abrupt stepwise jump of the caloric curve corresponding to the melting temperature of a 
bulk. 
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FIG. 2: Caloric curve for the pure Nii47 cluster. Temperatures T <T freezing and >T^melting 
correspond to the completely frozen and melted states, respectively. 
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FIG. 3: Heat capacity for the pure Nii47 cluster as a function of T. Temperature = 750 K 
indicates the melting temperature of the cluster (defined as the maximum of the heat capacity) . 
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The melting process can also be easily recognized by the peak in the temperature depen- 
dence of the heat capacity, as seen in Fig. [31 We have found that for the pure Nii47 cluster 
the melting temperature associated with the maximum of the heat capacity = 750 K 
is considerably smaller than the melting temperature for the bulk nickel T^""^ = 1728 K.— 
As was discussed above, the decrease in the melting temperature of the finite size clusters 
in comparison with the bulk occurs due to a substantial increase in the relative number 
of weakly bounded atoms on the cluster surface. According to the so-called 'Pawlow law' 
the melting temperature of a spherical particle possessing a homogeneous surface decreases 
linearly with increasing its inverse radius:- 



where R is the radius of the spherical particle and a is a constant, which can be defined 
by fitting the temperature Tm to available experimental data. Equation (jlj) is valid for the 
relatively large particles when shell effects are not important. The similar reduction in the 
melting temperature of the large Ni^v clusters with the number of atoms ranging from 336 

r 

to 8007 has been reported in Ref. 5. Based on the calculated melting temperatures from Ref. 
y, we have estimated the values of T^"''^ and a to be 1590 K and 3.65 A respectively. The 
radius of the cluster has been defined as R = VgN^^^, where = 1.375 A is the Wigner-Seitz 
radius for bulk nickel.— 

Apart from the main peak, the heat capacity displays an additional maximum at T ^ 
660 K — suggesting a stepwise melting process. This additional maximum corresponds to a 
slight change in the slope of the caloric curve beginning at T freezing ~ 600 K. Such a situation 
corresponds to the so-called pre-melted state — when the cluster surface melts first while the 
core of the cluster remains frozen. This is the stationary state where the coexistence of 
two phases — liquid surface and frozen cluster core — is observed. The exact delimitation of 
the two phases is relative but can be defined from the difference in mobility of atoms taken 
from the cluster surface and its core. Visualization of the molecular dynamics trajectories 
of atoms located at the surface and in the core of the cluster confirms this interpretation. 

Figures |4](a), |4](b), andlD^c) present the time dependence of the instantaneous values of 
the total energy calculated for the Nii47 cluster at the temperatures 600 K, 750 K, and 
800 K, respectively. The total simulation time is 10 ns. 

It is seen from Fig. |l](a) that the instantaneous values of the total energy oscillate around 
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FIG. 4: Time dependence of the instantaneous values of the total energy calculated for the Nii47 
cluster: (a) T = 600 K, frozen state; (b) T = 750 K, phase transition; (c) T = 800 K, molten state. 
The dashed-dotted lines demonstrate the time-averaged total energy of the system. 
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their time-averaged value at the cluster temperature T = 600 K. The fluctuations in the 
total energy of the system is a result of choosing the canonical ensemble of particles. To 
reproduce this ensemble, the temperature of the system is controlled by the Nose-Hoover 
thermostat, which produces the appropriate energy fluctuations.—'^ The frozen phase is 
characterized by thermal vibrations of the atoms around equilibrium positions in the cluster 
lattice, while the overall cluster structure remains unchanged. 

For the temperature T = 750 K, corresponding to the maximum of the heat capacity, 
the time dependence of the instantaneous values of the total energy changes, as shown in 
Fig. 111(b). The total energy of the system initially oscillates around its typical averaged 
value for the frozen Nii47 state. After some time, the system jumps to a disordered molten 
state — causing the total energy of the system to oscillate around the average value typical 
for the molten Nii47 state. The system then, after a finite period in the molten state, 
jumps back to its initial frozen state. This behavior is repeated and the system continously 
oscillated between the frozen-molten states. Note that, in the interval of temperatures 
T freezing < T < T melting, it is ucccssary to perform the molecular dynamics simulations for 
a relatively long time (on the order of 10 ns) to achieve time independence of the averaged 
total energy ( Etot ) of the system. The average life-time of the system in the frozen and the 
molten states depends upon the temperature of the cluster, resulting in a smooth change of 
the ( Efot ) from the frozen to the molten state as a function of temperature. Thus, Fig. IH^b) 
clearly demonstrates coexistance of the two thermodynamic phases of the finite system at 
the temperature of phase transition. A similar behavior of the time dependence of the short- 
time potential energy averages calculated for the Ari^ cluster at the temperature 33 K has 



been reported in Refs. 



72 



73 



Finally, for the temperature T = 800 K and above, the total energy of the Nii47 cluster 
oscillates around its typical averaged value for the molten state. 



B. C-doped Ni^y clusters 

We now focus our study on the thermodynamic properties of the C-doped Nii47 cluster. 
We consider a cluster with the C impurity in the vicinity of the cluster center as shown in 
Fig. [1Kb). 

Dashed-dotted lines in Figs. [5] and M demonstrate the temperature dependence of the 
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caloric curve and the heat capacity calculated for the C-doped Nii47 cluster. For the sake 
of comparison, we have plotted the same dependencies calculated for the pure Nii47 cluster 
in Figs. [5] and [6] by solid lines. As seen in Figs. [5] and [6l doping of the Nii47 cluster with 
a single C impurity reduces its melting temperature by 30 K. Thus, doping of a cluster 
consisting of more than a hundred of atoms by just one additional atom of impurity results 
in a considerable change of its melting temperature. 
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FIG. 5: Caloric curve for the pure Nii47 cluster (solid line) and the C-doped Nii47 cluster (dashed- 
dotted line). 

Figure \5\ demonstrates that the first change in the slope of the temperature-dependent 
( Etot ) calculated for the C-doped Nii47 cluster occurs at a temperature slightly below 600 K. 
This fact indicates that the C impurity does not have much influence on the surface melting 
of the cluster. Indeed, the C impurity is located near the cluster center where it induces 
local deformations. However, it does not significantly change the structure of the cluster 
surface, see Fig. dJ^b); therefore, it does not have much of an effect on surface melting. 

Further increase in the cluster temperature results in the melting of the C-doped Nii47 
cluster with the melting temperature = 730 K. The latter is defined by the position 
of the maximum of the temperature-dependent heat capacity. The decrease in the melting 
temperature of the C-doped Nii47 cluster can be explained as a result of the local distortions 
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FIG. 6: Heat capacity for the pure Nii47 cluster (solid line) and the C-doped Nii47 cluster (dashed- 
dotted line) as a function of T. Temperatures = 750 K and Tm = 720 K indicate the melting 
temperatures of the pure and the C-doped Nii47 cluster, respectively. 

of the cluster's icosahedral structure near the impurity. The C atom introduces a strain in 
the cluster lattice which decreases stability of the cluster and thus, its melting temperature. 

This effect is in accordance with the behavior of bulk material where it is known that 
defects in crystalline lattice provide nucleation sites for the liquid state (see, e.g., Ref. 



59 



and references therein) . However such an effect strongly depends upon the type of impurity. 



As was demonstrated in Ref. l28|, doping of the Agss cluster with Cu and Ni impurities in the 
center (where the bond length between Ag atoms and atoms of impurity are considerably 
smaller than the bond length between Ag atoms^^i^) results in an increase of the melting 
temperature of the system. Therefore substitution of the central Ag atom with an atom of 
Cu or Ni impurity induces the contraction of the cluster lattice in the vicinity of the center, 
which can partially release the strain of the icosahedral structure.— Thus the small impurity 
in the center allows the cluster to relax toward a configuration with better interatomic 
distances, thereby increasing the cluster stability.—*^ 

Our results demonstrate that doping of the C atom in the vicinity of the center of the 
Nii47 cluster induces an additional strain in the cluster structure and, as a result, the melting 



13 



temperature of the cluster decreases. To study in detail how this induced strain can change 
the thermodynamic behavior of the C-doped Nii47 cluster, we have performed calculations 
of the melting temperature of the cluster for the set of different parameters p and tq in the 
Morse potential ([3]) — in other words, modeling the variation in the Ni - C interaction. 
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FIG. 7: (Color online) Morse potential for the Ni - C interaction with different values of parameters 
p and tq. Curve 1, p = 3.295, ro = 1.763 A (optimal values); Curve 2, p = 1.648, tq = 0.882 A 
(reduced bonding); Curve 3, p = 4.943, ro = 2.645 A (enlarged bonding). The depth of the 
potential well is kept constant em = 2.431 eV. 

Curve 1 in Fig. [7] presents the dependence of the Morse potential V^'^''"{r) on the inter- 
atomic distance r caclulated for the optimal values of the parameters Em, P, and tq. Curve 
2 in Fig. [7] presents the potential V^^~'-''{r) when the bond constant Tq and parameter p 
are reduced by the factor 0.5 in comparison to their optimal values; while curve 3 in Fig. [7] 
presents the interaction potential V^'^~'-''{r) when the parameters ro and p are enlarged by 
the factor 1.5. Hence, the potentials presented in Fig.[7]by curves 2 and 3 model the impurity 
effect with the local compression and expansion of the cluster lattice respectively. 

To illustrate the effect of variation in the Ni - C interaction on cluster structure we 
present, in Figs. [H]^a)-(c), the histogram of the radial distribution of the number of Ni atoms 
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FIG. 8: The number of Ni atoms AA^ at distances between r and r + Ar from the center of mass of 
the C-doped Nii47 cluster calculated for different values of the parameters p and tq: (a) p = 1.648, 
ro = 0.882 A (reduced bonding); (b) p = 3.295, tq = 1.763 A (optimal values); (c) p = 4.943, 
ro = 2.645 A (enlarged bonding); (d) pure Nii47 cluster for comparison. The radial interval Ar is 
0.001 A. 

located at distances between r and r + Ar from the center of mass of the C-doped Nii47 
cluster (where the width of the radial interval is Ar = 0.001 A). The change in the cluster 
structure, due to varying the parameters of p and ro in the Morse potential, can thus be 
clearly seen. For comparison. Fig. M^d) represents the histogram of the radial distribution 
of the number of Ni atoms in the pure Nii47 cluster. 

In the case of the reduced bonding between the C impurity and the Ni atoms, the overall 
icosahedral shell structure of the cluster remains preserved — although there exists some 
relaxation of the lattice. Fig. [8]^a). The increase of the effective radius of the Morse potential 
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for the Ni - C interaction results in strong distortions of the icosahedral shell structure of the 
cluster, Figs.lHt^b) and[8](c). Such distortions reduce the stability of the cluster in comparison 
to the compact icosahedral structure. 
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FIG. 9: (Color online) Caloric curves for the C-doped Nii47 cluster calculated for the different 
values of parameters p and rg. Curve 1, p = 3.295, rg = 1.763 A (optimal values); Curve 2, 
p = 1.648, ro = 0.882 A (reduced bonding); Curve 3, p = 4.943, tq = 2.645 A (enlarged bonding). 
The depth of the potential well is kept constant em = 2.431 eV. 

Figures M and dU] demonstrate the temperature dependence of the caloric curve and the 
heat capacity of the C-doped Nii47 cluster calculated for the different parameters of the 
Morse potential V'^^~'"{r). Curves 1 in Figs. [9] and [TO] present the caloric curve and the 
heat capacity calculated with the optimal values of parameters for the Ni - C interaction. 

Decreasing the effective radius of the Morse potential (see curve 2, Fig. [7]) by a factor 0.5 
(c./. the optimal values) results in the contraction of the cluster lattice in the vicinity of 
impurity. This contraction releases the strain in the icosahedral structure, therefore leading 
to an increase of the cluster stability. In this case the, melting temperature of the impurity 
doped Nii47 cluster increases by 21 K, compared the pure cluster (see curve 2 in Fig. [TOj) . 

An increase in the effective radius of the Morse potential for Ni - C interaction to its 
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FIG. 10: (Color online) Heat capacity for the C-doped Nii47 cluster calculated for the different 
values of parameters p and tq. Curve 1, p = 3.295, ro = 1.763 A (optimal values); Curve 2, 
p = 1.648, ro = 0.882 A (reduced bonding); Curve 3, p = 4.943, ro = 2.645 A (enlarged bonding). 
The depth of the potential well is kept constant em = 2.431 eV. Arrow indicates the melting 
temperature of the pure Nii47 cluster. 

optimal value (see curve 1 in Fig. [7]) results in the appearance of an additional strain in the 
cluster lattice. Therefore, the melting temperature of the C-doped Nii47 cluster decreases 
compared to the pure Nii47 cluster. 

The further increase in the effective radius of the Ni - C interaction creates a strong 
deformation and rearrangement in the cluster structure, see Fig. [8](c). In this case, location 
of the impurity in the vicinity of the cluster center becomes energetically unfavorable. Thus, 
the impurity atom shifts towards the cluster surface to minimize the destruction of the 
icosahedral lattice. This effect prevents the further decrease in the melting temperature of 
the cluster, as an impurity located on the cluster surface will not have much influence on 
the melting transition of the core. 

Melting transition can be also recognized from the analysis of the trajectories of the 
individual atoms and their diffusion in the volume of the cluster. The melting transition 
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occurs when atoms begin their Brownian motion instead of the thermal vibrations around 
their equihbrium positions in the ordered cluster's lattice. Such a transition can be seen as 
a step in the temperature dependence of the diffusion coefficient. The diffusion coefficient 
of the atom i in a media is defined as (see, e.g., Refs. 00): 

a4^<.?W>. (5) 
where {rfit)) is the mean-square displacement averaged along the atom's trajectory 

^ nt 

W> = - E (I^^(^o^- + ^) - ^^^(^0,))' . (6) 

Here Ri()f:) is a radius vector of an atom i at the time t, and rit is the number of time origins, 
toj, considered along the trajectory. Equation ([5]) is valid for the simulation time smaller 
than the time required for a particle to migrate across the diameter of the cluster. 

Figure [TT] demonstrates the 2D projection of trajectories calculated for the Ni atom in the 
center of the C-doped Nii47 cluster (filled dots), the Ni atom from the vertex of the cluster 
surface (stars) and the C impurity (open dots). The temperature of the cluster ranges from 
400 K to 800 K, as is shown in Fig. [TTl The output time step is 1 ps with the total simulation 
time of 2 ns. 

As seen in Fig. [TT|, for a low temperature, T = 400 K, all selected atoms in the C-doped 
Nii47 cluster vibrate around their equilibrium positions. At this temperature, the cluster 
is frozen — as per the analysis of the caloric curve and the heat capacity calculated for the 
C-doped Nii47 cluster, see Figs. |5]and[6]). 

It has been discussed above that the temperature dependence of the heat capacity of the 
C-doped Nii47 cluster exhibits two maxima corresponding to surface and volume melting of 
the cluster. It is seen from Fig. E] that the first maximum in the temperature dependence of 
the heat capacity appears when the cluster is at T= 640 K (surface melting), while the second 
maximum appears at T= 720 K (volume melting). The temperature of 680 K corresponds 
to the intermediate state when the surface of the cluster has already melted but the cluster 
core is still frozen, as confirmed by the analysis of the atomic trajectories. As can be seen in 
Fig. [TT]at T= 680 K, the surface Ni atom begins to diffuse on the surface, while the central 
Ni atom and the C impurity are still vibrating around their equilibrium positions. Note, 
that the icosahedral surface is inhomogeneous and consists of 12 vertices, 20 faces and 30 
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FIG. 11: (Color online) The 2D projection of trajectories calculated for the Ni atom in the center 
of the C-doped Nii47 cluster (filled dots), the Ni atom from the vertex of the cluster surface (stars) 
and the C impurity (open dots) for the cluster temperatures T = 400 K, 680 K, and 800 K. The 
output time step is 1 ps and the total simulation time is 2 ns. 
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edges. The binding energies of the atoms taken from the vertices, faces and edges are thus 
shghtly different. Hence, these atoms begin diffusion at different temperatures. 

Finally at the temperature of 800 K, cluster has completely melted as can be observed 
from the temperature behavior of the heat capacity presented in Fig. El Figure [H] demon- 
strates that the Ni atoms are moving in the entire volume of the cluster. The C impurity 
is also moving randomly in the cluster volume although the movement only occurs in the 
central part of the cluster. Hence, this supposes the heterogeneous distribution of the C 
atoms in the melted Nii47 clusters at T= 800 K. 

The change in cluster structure upon melting is clearly seen in the time-averaged radial 
distribution of atoms in the cluster. The radial atomic distribution function g{r) is defined 
as: 

g{r) = dN{r)/dr, (7) 

where dN{r) is the number of atoms in the spherical layer at distances between r and r + dr 
from the center of mass of the cluster. 

Figure [T2] demonstrates the time- averaged radial distribution function g{r) for Ni and C 
atoms in the C-doped Nii47 cluster calculated at cluster temperatures of T = 300 K, 600 K, 
700 K, 720 K, 800 K, and 1000 K. The chosen range of temperatures allows for the analysis 
of the cluster structure in the frozen, transitional and molten states. 

At the low cluster temperature, T = 300 K, one can see the icosahedral shell structure of 
of the nickel subsystem of the C-doped Nii47 cluster, consisting of the central atom and three 
icosahedral shells. The second and third shells are split — corresponding to atoms in vertex 
(12 atoms per shell) and non- vertex positions. Heating the cluster up to 600 K washes out 
the subshell splitting, nonetheless the icosahedral shells remain well separated. 

With increase in the cluster temperature, up to 700 K, the second and the third shells 
beging to merge, although the first and the second shells are still separated (i.e., the radial 
distribution function g{r) is equal to zero in the space between shells). At the temperature 
corresponding to the maximum in the heat capacity of the C-doped Nii47 cluster, T = 720 K, 
the first and second icosahedral shells merge. 

Finally, at the temperatures corresponding to the molten state (800 K and 1000 K in 
Fig. fT2l) . the distribution of Ni atoms become more homogeneous and the sharp shell struc- 
ture washes out. Nevertheless, even at T = 1000 K, some radial order with maxima at 2.6 
A, 4.6 A, and 6.5 A still remains, suggesting that even a molten cluster of a finite size 



20 






800 K 




C atom 


1000 K 

























3 



FIG. 12: (Color online) The time-averaged radial distribution function g{r) calculated for the Ni 
and C atoms in the C-doped Nii47 cluster at different temperatures. 
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manifests some signs of a shell structure. This effect might be a general feature of finite 
systems, similar to the surface-induced ordering in liquid crystals or the layering effect at 
free liquid surfaces.— "^"^ 

Figure fT2l demonstrates that, at cluster temperatures 300 K and 600 K, the time-averaged 
radial distribution function calculated for the C impurity atom exhibits a sharp maximum 
at distances 3 A from the center of mass of the C-doped Nii47 cluster. Thus, at these 
temperatures, the C impurity is located between the first and the second icosahedral shells 
of the Ni atoms. By further increasing the cluster temperature to temperatures near the 
phase transition region, the radial distribution g{r) of the C impurity becomes wider and 
the appearance of a second maximum at distances ^ 2 A can be observed. While at T > 
800 K, the C impurity can be found to be distributed in the central part of the cluster. 
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FIG. 13: (Color online) Temperature dependence of diffusion coefficients calculated for the Ni atom 
in the center of the C-doped Nii47 cluster (filled dots), the Ni atom on the cluster surface (stars) 
and the C impurity (open dots). 

Figure [13] demonstrates the temperature dependence of diffusion coefficients calculated 
for the selected atoms in the C-doped Nii47 cluster. It is seen from Fig. [13] that the Ni atom 
on the cluster surface begins to diffuse at temperature of 600 K, while Ni atom in the cluster 
center and the C impurity remain frozen. 
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As has been shown above, the caloric curve calculated for the C-doped Nii47 cluster 
(see Fig. [5]) clearly demonstrates a stepwise melting behavior. In the temperature interval 
600 < T < 700 K, the slope of the caloric curve changes slightly in comparison with that 
for the frozen state, suggesting the existence of a pre-melted state. For the temperatures 
700 < T < 760 the time-averaged total energy ( Etot ) growing rapidly and finally reaches 
its asymptotic behavior at T > 760 K. The variations in the slope of the caloric curve 
result in the appearance of the two maxima in the temperature dependence of the heat 
capacity. These maxima are associated with the surface and core (volume) melting of the 
C-doped Nii47 cluster. This proposition is fully confirmed by the analysis of the temperature 
dependence of the diffusion coefficients for the Ni atoms from: a) the cluster surface and b) 
the cluster center. 

Thus, in the temperature interval of 600 < T < 700 K, only surface atoms diffuse, 
confirming that the initial melting of the cluster surface. The Ni atom located in the cluster 
center then begins to diffuse at T= 700 K. The diffusion coefficients calculated for Ni atoms 
located on the cluster surface and in the core are different up to T = 760 K. The difference 
in the diffusion coefficients show that the surface and the core atoms are not yet fully mixed 
in the cluster volume. The difference disappears at cluster temperatures of T > 760 K, when 
the C-doped Nii47 cluster has become completely molten. 

Figure [13] demonstrates that the C impurity begins to diffuse at T ^ 700 K — similar 
to the central Ni atom. However, values of the diffusion coefficients calculated for the C 
impurity are considerably lower than those for the Ni atoms. Knowledge of the diffusion 
coefficients of the C impurity in a nickel cluster can be used for building a reliable kinetic 
model of carbon nanotube growth.—*^ 

IV. CONCLUSION 

Doping of Nii47 with a carbon impurity lowers its melting temperature by 30 K due to 
excessive stress on the cluster lattice. The magnitude of the change induced is dependent 
upon the parameters of the interaction between the nickel atoms and the carbon impurity. 
We have demonstrated that an induced contraction of the icosahedral cluster's lattice in 
the vicinity of the impurity results in an increase of the melting temperature of the cluster; 
whereas additional strain in the lattice results in the reduction of the melting tempera- 
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ture. Therefore, the melting temperature of atomic clusters can be effectively tuned by the 
addition of an impurity of a particular type. 

Doping by a C impurity changes the melting temperature of the cluster, consequently 
this means that doping affects the mobility of the atoms in the Ni cluster. This effect has 
to be taken into consideration in particular applications with metal clusters when the entire 
process depends on the thermodynamic state of the cluster. An example of such experiment 
is the process of the catalytically activated growth of carbon nanotubes. The kinetics of 
the carbon nanotube growth depends upon diffusion of carbon atoms through the metal 
catalyst. Presence of the impurities can considerably change the flux, thereby affecting the 
growth rate of the carbon nanotube. The additional change in the thermodynamic state of 
the catalytic particle in the nanotube growth process might also depend on the strength of 
the interaction of the particle with a substrate. 

In the present work, we have considered a single C impurity in the cluster of Nii47. We in- 
tend to further study the effect of how several C impurities will influence the thermodynamic 
properties of the host cluster. In particular, it is important to find the optimum conditions 
(concentration of C atoms, temperature, thermodynamic state of the particle, etc.) when 
the C atoms begin aggregating into ordered carbon structures, such as nanotubes. 

The influence of impurities on properties of finite systems is a general effect. While 
our results were obtained for free clusters, many interesting problems can be found when 
one considers the influence of impurities on the phase transitions and stability of clusters 
deposited on a substrate. Thus, recently it has been experimentally shown that the oxidation 
of silver clusters deposited on a HOPG surface changes the stability and morphology of 
cluster formations.—"^ Clusters on substrates have important technological applications and 
the understanding of how these clusters stabilize on the surface are of profound interest. 
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